Realizability of stationary spherically symmetric transonic accretion 
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The spherically symmetric stationary transonic (Bondi) flow is considered a classic example of an 
accretion flow. This flow, however, is along a separatrix, which is usually not physically realizable. 
We demonstrate, using a pedagogical example, that it is the dynamics which selects the transonic 
flow. 
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In astrophysics, the importance of accretion processes 
can hardly be overstated, especially in the context of the 
study of compact astrophysical objects and active galac- 
tic nuclei. Of such processes, a very important paradigm 
is that of the steady spherically symmetric flow, in which 
the motion of the accreting matter is steady and spher- 
ically symmetric, obeying the boundary condition that 
the bulk velocity is to fall off to zero at infinity, while the 
density asymptotically approaches a fixed value. Since 
this is the simplest situation in the class of accretion 
flows, it is the starting point in all relevant texts 
Studied extensively, among others, by Bondi about 
fifty years ago, it also gives one of the clearest examples 
of a transonic fiow (with which, in spherical symmetry, 
Bondi's name is associated), i.e. a flow in which the ve- 
locity is subsonic far away from the star and becomes 
supersonic as the surface of the star is approached. It is 
this simplicity and easy comprehensibility of the model, 
that overshadows the fact that it is not encountered most 
often in practice. Accordingly, the sphericallysymmetric 
flow has been of continuing interest |j, Isl lo, LD| ■ We first 
argue that the transonic (Bondi) solution of the steady 
spherically symmetric flow, would not be physically re- 
alizable. This proposition is in contradiction to the po- 
sition of Bondi himself, that "the case physically most 
likely to occur is that with the maximum rate of accre- 
tion" 3], which, in spherically symmetric accretion, is a 
case that is readily identified as the transonic solution 
0. We then argue that it is the dynamics which actu- 
ally selects the transonic flow from among all possible 
trajectories. 

The typical solutions for the spherically symmetric flow 
in the velocity-coordinate space are shown in Fig^ (to be 
read without the arrows). The two dark solid curves la- 
belled A and W refer to the accretion flow and the wind 
flow respectively. The intersection point is an equilib- 
rium point. The problem with Figgis that when it is 
read without the arrows it is slightly misleading. It does 
not show along what route an integration of dv/dr would 
proceed if we start with an initial condition v — Vm at 
T = Tin far away from the star. For a physically realiz- 
able flow, an initial condition inflnitesimally close to a 



point on the accretion curve A, would trace out a curve 
inflnitesimally close to A and in the limit would correctly 
reproduce A, evolving along it and passing through the 
equilibrium point (sonic point) as we integrate dv/dr, 
obtained from Euler's equation. We will show that the 
arrows on the integration route are as shown in Fig^ 
The direction of the arrows indicates that the spherically 
symmetric transonic accretion flow should not be physi- 
cally realizable. 

Is this a result confined to the spherically symmetric 
flow? We have verified, as we shall show briefly, that 
this holds also for the axisymmetric rotating accretion 
flow for a thin disc which is a situation of practical 
interest. The technique for assigning the arrows remains 
identical to that of the spherically symmetric case . 

We wish to make our point with a fairly straightfor- 
ward example. We consider the differential equation 



y 



dy _ Y{x,y) _ 

dx X{x, y) y — X 



(1) 




Fig.l 



FIG. 1: Phase trajectories for spherically symmetric accre- 
tion onto a star. The bold solid lines, A and W, represent 

"accretion" and "wind" respectively. The fixed point is at 
,.2 , 



r = ro and v /cs 



1. Linear stability analysis indicates 



that the fixed point of the flow is a saddle point. The direc- 
tion of the arrows along the line A, demonstrates that the 
transonic flow is not physically realizable. 
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Fig. 2 



FIG. 2: Integration of Eq.Q gives a pair of straight lines, 
with the integration constant fixed by the intersection point 
(1,1). In the figure the lines are marked A' and W'. 
Parametrizing Eq.Q as Eq.Q, we get the concept of arrows. 
Linear stability analysis indicates that for A' and W', the 
intersection point (1, 1) is actually a saddle point, for which 
the arrows are as shown above. 



whose integral can be written down as 



- 4x + 2xy = -C 



(2) 



where C is a constant. If we want that particular so- 
lution which passes through the point where Y{x,y) = 
X{x,y) = 0, namely x = y = 1, then C = 2. The curve 
x"^ — y"^ — Ax + 2xy = —2 factorizes into a pair of straight 
lines : ?/-x(y2+l) + V2 = and y + a;(V2-l)-\/2 0. 
This pair is shown in Fig|21 (to be read again without the 
arrows) as the lines marked A' and W. 

We want to explore the process of drawing the line 
A' from a given initial condition. On line A', y = at 
X = 2 + \/2. Let us begin with the initial condition y = 
at X — 2 + ^/2 — €, where < e ^ 1 . The given initial 
condition fixes the constant C as C = 2[1 + \/2e — e^/2]. 
Using this value of C, we can plot the curve given by 
Eq.(|2Jl. For a given x, a value of y is given by the relevant 
root of the quadratic equation thus obtained. The two 
roots are 



■±V2 



ix-l)' + V2e-- 



-,1/2 



(3) 



from which it is clear that to satisfy y = 0a,tx = 2 + 
V2 — e, the negative sign has to be chosen to let us have 



y^x 



V2 



{x - if + V2e 



1/2 



(4) 



an error of an infinitesimal amount e in prescribing the 
initial condition on A', i.e. if we prescribe y = at 
X = 2 + \f2 — e instead of?/ = 0atx = 2 + V2, then 
the "error" made at a; = relative to A' is 2\/2 which is 
0(1). An infinitesimal separation at one point leads to a 
finite separation at a point a short distance away. This is 
what we mean by saying that line A' (and similarly W) 
should not be physically realized. 

The clearest and most direct understanding of the 
difficulty is achieved by writing Eq.lQ} as a set of two 
parametrized differential equations 



dy_ 

dT 
dx 

d7 



x + y-2 



y-x 



(5) 



where r is some convenient parameter. The fixed point of 
this dynamical system is (1, 1), namely the point where 
Y{x, y) and X{x, y) are simultaneously zero — the point 
through which both A' and W pass. Linear stability 
analysis of this fixed point in r space shows that it is a 
saddle. The critical solutions in this y — x space can 
now be drawn with arrows and the result is as shown in 
FigEl The distribution of the arrows (characteristic of 
a saddle implies A' and W cannot be physically 
realized. 

For the spherically symmetric flow, the relevant vari- 
ables are the radial velocity v and the local density p. 
Ignoring viscosity, we write down Euler's equation for v 
as 



dv 



dv 



dr 



1 dP 
p dr 



dV_ 

dr 



(6) 



where P is the local pressure and V = V{r) is the po- 
tential due to a gravitational attractor of mass M, i.e. 
V{r) ~ —GM/r. The pressure is related to the local 
density through the equation of state P — Kp"' where K 
is a constant, and 7 is the polytropic exponent with an 
admissible range given byl<7<5/3. The local density 
evolves according to the equation of continuity, 



l^ + l2^(/'^^ J 



dt 







(7) 



The stationary solution imphes dv/dt = dp/dt = 0. The 
local sound speed is given by Cg^ = dP/dp = ^Kp^^^ 
and in the stationary situation, we can use Eq.((7J) to 
write Eq.® as 



At a; = 0, we then get y = -\/2(l + ^/2e - eV2)^/^ 
which is very different from y = y/2, that one gets on line 
A'. In the limit of e — > 0, one generates a part of A' 
(x ^ 1) and a part of W {x ^ 1), instead of the entire 
line A'. Another way of stating this is the sensitivity 
to initial conditions in the drawing of A'. If we make 



d{v^) 2w2 (2c,2 - GM/r) 



dr 



(8) 



whose integration with different initial conditions is sup- 
posed to generate the curves in Fig^ To assign arrows 
we write Eq.® in a parametrized form 



3 



dr 
dr 
d7 



2v^ 2c. 



GM 



(9) 



to find that the critical point is at r = rp and v — vq, such 
that vq^ = Cso^ and Cso^ = GM/2rQ. It is immediately 
clear that the critical point is the so called sonic point, 
and it can be fixed in terms of the constants of the sys- 
tem with the help of suitable boundary conditions . We 
now carry out a linear stability analysis around the fixed 
point by writing = VQ^il + Si) and r = ?'o(l + (52)- Lin- 
earizing m 5i and 82, followed by some straightforward 
algebra then yields 



d5i 

dh 
dr 



CsO 



- (7 - l)6i + (6 - 47)<52 



(10) 



The solutions for 5i^2 are of the form e'*''^, where A is 
to be found from the roots of 



A-|-2(7-l)cso2 2(47-6)cso2 
-[(7 + l)/2]c,o^ A-2(7-l)c,o2 



= 



(11) 



The eigenvalues A are given by A = ±Cso^-\/2(5 — 87). 
For 7 < 5/3, the fixed point is a saddle. Since the physi- 
cal situations are restricted to this range, we find that the 
fixed point is always a saddle and hence the distribution 
of arrows will be as in Fig^ 

For the more realistic axisymmetric case (as in a thin 
accretion disc), using R as the radial distance in the 
plane, the static equations of the flow for an angular ve- 
locity Q. and for the disc thickness , are 

The Equation of Continuity : 



pRHv ~ constant (12) 
The Momentum Balance Equation : 



1 div^) 



2 dR pdR 
The Angular Momentum Balance Equation : 



'"dR^ ^ " pRH dR \ dR 



(13) 



(14) 



where = GM/R^ and a is the effective viscosity of 
Shakura and Sunyaev j^. An integral of motion follows 
from Eg. 1(1^ and we can obtain an equation for d{v'^)/dR 
akin to Eq. , which we can write as the dynamical sys- 
tem 



d{v^) 
dr 
dR 
d7 



2v' 



R 



1 + 7 



1-^7 



(15) 



from which, taking the inviscid limit (i.e. a = 0), an 
analysis identical to the one following Eq. (O leads to the 
eigenvalues A of the stability matrix. Of these one can 
easily be shown to be appropriate for a saddle 0] . 

We now focus on the fact that the real physical problem 
is dynamic in nature as exhibited in Eqs.® and lO, even 
as we have tacitly assumed that we can study the prob- 
lem directly in its stationary limit. As it turns out there 
are an infinite number of stationary solutions and if the 
static limit is taken directly, the Bondi solution becomes 
very sensitive to the process of integration to determine 
the stationary trajectory. Several studies in the linear 
stability analysis of the stationary flows that have been 
carried out in the past have not really helped clarify the 
situation as it was finally established that among the sta- 
tionary solutions, the transonic Bondi solution as well as 
all the subsonic solutions are stable in the linear stability 
sense 0, ^3 • -^^Y selection mechanism then must be of 
a non-perturbative nature. 

Accordingly, we return to our pedagogic example of 
Eq.(^ but now consider ?/ as a field y{x, t) with the evo- 
lution 



^ + {y-x)-=y + x-2 



(16) 



The stationary solutions y{x) satisfy Eq.Q and Fig[21 
shows that their separatrices are given by y — a; (a/2 + 
1) + = and 2/ -t- x{V2 - l) - %/2 = 0. It should 
be easy to see that a linear stability analysis in t around 
the family of stationary solutions y{x) would show an 
infinite number of them to be stable, and therefore each 
one of them could be a perfectly valid physical solution. 
We will now show that the non-perturbative dynamics 
actually selects the separatrices. 

The general solution of Eq.ll(j|l can be obtained by 
the method of characteristics 12]. The two independent 
characteristic curves of Eq. H16() are obtained from 



dt 
T 



dx 



dy 



y — x y + X ■ 



(17) 



They are 



4 



- 2xy -x^ +Ax^C 



c 



(18) 



in which the latter is obtained by integrating the dx/dt 
integral with the help of the first solution. With the 
choice of the upper sign now (since we should want the 
evolution to proceed through a positive range of i), the 
solutions of Eq. H16(l can be written as 



-2xy-x'^ +Ax = / 



x/2 



'V2t 



(19) 



where / is an arbitrary function, whose form is to be 
determined from initial conditions. In this case we choose 
the initial condition that y{x) = at i = for all x. This 



leads to f{x[l- 1/V2] - 1) = -x^ 
for / as f lz) = -2{2V2 + 3)^2-4 
solution to Ea. (|16|) is then given as 



Ax to give a form 
2 + l)z + 2. The 



y-x{V2 + l) +V2 y + x{V2-l)-V2 
= -2(2 + V2) \x{V2 -l)+y-V2 
-(3 + 2%/2) [a;(V2- 1) + y - V2 



(20) 



Clearly for t — > cx), the right hand side in Eq. 12()|l 
tends to zero and we approach one of the two separatri- 
ces shown in FigEl Of the two possible separatrices, the 
one which will be relevant will be determined by some 
other requirement. For the astrophysical flow the two 
separatrices are the Bondi accretion flow and the tran- 
sonic wind solution. One chooses the proper sign of the 
velocity to get the flow in which one is interested. 

The mechanism for the selection of the asymptotes in 
FigH ^ the favoured trajectories is identical. This can 
be appreciated from a look at Eq.|(Bl). In the "pressure- 
free" approximation, we have a set of stationary solutions 

— 2GM/r — (? . Which of these would be selected by 
the dynamics? If we start from u = at t = for all r, 
an identical reasoning to the one given above shows that 
it is the path with c = which is selected. The solution 
of the differential equation 



dv 
9t 



dv 



GM 



dr 



by the method of characteristics yields 



GM 



GM (v 



exp 



' 2rv 



(21) 



^) (22) 



in which = 2GM/c^ . Obviously for t — > oo we will 
get the static solution = 2GM/r. This is also the 
choice according to the energy criterion. Corresponding 
to the initial condition given, the lowest possible total 
energy is = /2 — GM/r ~ and the dynamics se- 
lects this particular stationary trajectory. This selection 
mechanism is entirely non-perturbative in character and 
provides the mathematical justification for Bondi's asser- 
tion that the energy criterion should select the stationary 
flow With the pressure term included, this argument 
leads to the selection of the transonic path i.e. gravity 
always wins over pressure at small radial distances. 
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